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\Q ■ Abstract 



A method is suggested to deduce the anisotropy in neutral pions 



■ by measuring the azimuthal anisotropy of photons in ultra-relativistic 

CD nuclear collisions. The ratio of the estimated anisotropy in photons to 

the anisotropy in neutral pions is seen to scale with a parameter which 



depends on photon multiplicity and anisotropy. This parameter can be 
determined from experimental data. 



1 Introduction 

In ultra-relativistic nuclear collisions the colliding nuclear matter emerges with 
a pattern that has its origin in the incompressibility of nuclear matter and the 
mean field effects. The observation of collective behaviour in the system vali- 
dates the use of a hydrodynamic description of the evolution of the collision. 
Within the hydrodynamic model the collective motion results from pressure 
gradient in the matter, which depends upon the compressibility of the under- 
lying equation of state [|IJ. In the case of a phase transition from hadronic 
matter to quark gluon plasma, increased number of degrees of freedom may 
lead to a softening of the equation of state @, p|, U, ||. Collective flow mani- 
fests in altered event shapes which can be studied by measuring the momentum 
distributions of produced particles. The azimuthal distribution of these event 
shapes are studied using coefficients in their Fourier expansion. One of the 
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major current aims of the study of ultra-relativistic heavy ion experiments is 
to determine the magnitude of anisotropy present in the final state momen- 
tum distributions of various particle species, in different regions of phase space. 
The systematic dependence of this anisotropy on energy, centrality and other 
collision parameters is predicted to be sensitive to the equation of state [0. 

Most of the experimental studies on azimuthal anisotropy have been re- 
stricted to charged particles (a large number of results for charged particle 
measurements can be found in the review article |§). The measured anisotropy 
of charged particles includes effect due to final state Coulomb interactions and 
may inhibit a direct comparison of results with the predictions of strong inter- 
action physics. 

Measurements of anisotropy in the azimuthal distribution of photons have 
been reported recently || 0]. Since almost 90% of photons produced in ultra- 
relativistic nuclear collisions originate from the decay of 7r°, the anisotropy 
measured in photons should have their origin in the anisotropy of it produced 
in the reaction. However, the measured anisotropy in photons will also in- 
clude some effect due to decay. The decay introduces non-flow correlations 
amongst kinematic variables of photons due to four-momentum conservation 
and also dilutes the existing correlation in azimuthal angles. Determination 
of the effect of decay should enable deduction of the anisotropy in n°. This 
would complement the study of flow using charged particles and would be 
free from Coulomb effects. Such a determination will provide much needed 
additional information to characterise the final state momentum distributions 
[§]. A comparison of the anisotropies deduced from charged and neutral pions 
would help in estimating the effect of final state Coulomb interaction on the 
charged particles. 

In the present work we study the effect of 7r° decay on the observed 
anisotropy of photons for various assumed flow parameters and multiplicities 
using a fast simulation technique. Preliminary results can be found in || ||], |10[ . 
We show how the anisotropy parameter for the parent 7r° distribution can 
be deduced using the measured anisotropy in the azimuthal distribution of 
photons. 



2 Simulation 



A complete determination of the flow parameters in different regions would 
require a measurement of the three-dimensional momentum distribution of all 
the particles. Several experiments determine the anisotropy parameters in az- 
imuthal angle distributions (dN/d(p) of produced particles. [|, 0, [Tl|, |12T|. In the 



present work we restrict ourselves to the study of photon multiplicity in a de- 
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tector having full azimuthal acceptance in one unit of pseudo-rapidity suitably 
chosen with practical detectors in mind J7|, We generate events having 

azimuthally anisotropic shapes and analyse them using methods employed for 
the analysis of experimental data. 



2.1 Anisotropy 

Let us assume that the event shape in the azimuthal plane can be described 
by an arbitrary ellipse. The most general equation of an ellipse in polar coor- 
dinates is 

r(0) = [1 + 2«i cos(0 - V>i) + 2v 2 cos 2(0 - V> 2 )] (1) 

where ipi and ip2 define the direction of the shift of the centroid of the ellipse 
from the origin and the orientation of the major axis of the ellipse respectively. 
v\ is the magnitude of the shift of the centroid and v 2 is a measure of the 
difference in the major and the minor axes of the ellipse. v\ and v% denote 
the first order and the second order flow components, termed as directed and 
elliptic flow. 



2.2 Event generation 

For the fast simulation, a set of events consisting of 7r + ,7r~ and 7r° is gener- 
ated. The particle multiplicity in the event is selected in such a way that the 
multiplicities of charged particles 7r~) and decay photons within the one 
unit of pseudo-rapidity acceptance of the given detector are in the range of 100 
to 4000. Such a wide range of multiplicity corresponds to varying centralities 
encompassing the SPS, RHIC and LHC energies and serves to determine the 
effect of fluctuations due to finite number of particles. The multiplicity distri- 
bution of particles is assumed to be Gaussian with a width a = 2\[M for a 
mean multiplicity M. 

We assume the azimuthal distributions to have initial flow components 
v™, where m (=1,2) denotes the order of flow. Using these values, the az- 
imuthal angles of the charged and neutral pions are assigned according to the 
distribution given by eqn. (|l|). ipi is chosen randomly. For the present work 
we have chosen ip2 = which corresponds to the major axis of the ellipse 
along the direction of the shift of the centroid. Physically this represents in- 
plane flow. The data are also generated for different values of vf 1 and v™, 
corresponding to different initial flow. 

The pseudo-rapidity and p^-distributions are taken from a parametrised 



form of the HIJING event generator |T3[ at the relevant energy. Neutral pions 
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are generated in the pseudo-rapidity interval extending to 1.5 units on both 
sides of the assumed acceptance of the detector. These are allowed to decay 
and the T], pt and of the photons are obtained. The photons falling within 
the acceptance region of the detector are used for further analysis. 



2.3 Analysis 

The set of events thus generated is analysed to extract the anisotropy param- 
eters v m using the Fourier expansion technique in a manner identical to the 
analysis of experimental data. The present analysis is based on 10000 events 
for each combination of multiplicity and anisotropy. 

The generated data have been analysed to obtain the following: 



1. ipl st : the estimated direction of the first order event plane. Because of 
finite particle multiplicity the estimated event plane i[)f st in the generated 
data fluctuates about the actual event plane, ipi, with a spread that 
depends on the initial anisotropy and the multiplicity. 

2. v[ = (cos(0 — ipl st )) '■ a measure of the directed flow in the direction of 
ipl st , where the average is over all particles of all events. 

3. ip2 St '■ the estimated direction of second order event plane. This also 
fluctuates about the actual direction ip 2 . 

4. v' 2 = (cos 2(0 — ip2 St )) '■ a measure of the ellipticity about ^f 5 * where the 
average is over all particles of all the events. This is a measure of the 
difference of the major axis and the minor axis of the event shape. 



Since the estimated event plane differs from the actual event plane, the 
values of v' m are systematically lower by a factor (cos m{jp m — ^')), where 
the average is over all events for a particular sample. This factor is the res- 
olution correction factor (RCF). The values v m that represent the magnitude 
of directed (m =1) and elliptic (m = 2) flow, for the sample of events, are 
obtained by v' m /RCF. 

Experimentally, RCF is obtained using the subevent method described 
in ||14|| . Here every event is divided randomly into two subevents of equal 
multiplicity and the angle ip m is determined for each subevent. This enables 
a determination of a parameter Xm directly from the experimental data using 



the relation |14|, H| : 



N events (m\r m -r m \>7r/2) 



ml 

Ntotal 



(2) 
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where N tota i denotes the total number of events, ip^ ty h m are the estimated 
angles of the two subevents (labelled a and b) and the numerator on the left 
denotes the number of events having the angle between subevents greater than 
n /2m. The parameter Xm is then used to determine RCF according to the 



relation given in the Ref. flTJ 



The generated azimuthal distribution of the charged pions are also anal- 
ysed in a manner similar to the azimuthal distribution of photons. Comparing 
the results of charged particles and photons enables a determination of the 
effect of decay. 




4000 2000 
Multiplicity 



Figure 1: Estimated value V\ of directed flow for charged particles and pho- 
tons as a function of multiplicity. Symbols denote different initial anisotropics 
varying from to 0.1. For photons the initial anisotropy was introduced in the 
7T° distribution. The error bars are statistical. 



3 Results 



3.1 Effect of decay 

The results for the estimated anisotropy v\ for directed flow are shown in 
Fig. |I](a) as a function of multiplicity for charged particles and Fig. |l|(b) for 
photons for varying values of v™. Similar results for elliptic flow are shown 
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in Fig. |2|. For charged particles the analysis technique correctly reproduces 
the initial anisotropy for all multiplicities for both directed and elliptic flow, 
demonstrating the validity of the technique used. For lower multiplicities the 
errors on the estimated values are larger due to greater fluctuations. 
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Figure 2: Estimated value t>2 of elliptic flow for charged particles and pho- 
tons as a function of multiplicity. Symbols denote different initial anisotropics 
varying from to 0.1. For photons the initial anisotropy was introduced in the 



7T 



-distribution. The error bars are statistical. 



The same analysis technique, when applied to photons, yields values of 
v m that are generally different from the initial anisotropy present in the 7r° 
distribution. This difference in the anisotropy values for photons and 7r° is 
attributed to two effects of decay : (i) it dilutes the initial flow by randomising 
the existing azimuthal correlation, and (ii) it mimics flow, the correlations 
due to momentum conservation producing an anisotropy in the final state 
distributions. 

To estimate the effect of momentum conservation correlation we gener- 
ated azimuthally symmetric distributions, corresponding to v™=0. The decay 
photons falling within the acceptance of the detector were analysed, and the 
values of v m obtained for them are also shown in Fig. |l](b) and in Fig. 0(b). One 
observes that the deduced anisotropy in photons is always non-zero. At low 
multiplicities the values can be as high as those observed in charged particle 
measurements flT2"|. The values for elliptic flow are systematically lower than 
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for directed flow, implying lesser effect of momentum conservation correlation 
on the elliptic event shape. 

For all values of initial anisotropy, the estimated anisotropy in photons 
decreases with increasing multiplicity. At low multiplicity the correlation due 
to momentum conservation competes with dilution due to decay. The rela- 
tive importance depends on the initial anisotropy in neutral pions. For low 
anisotropy in 7r°, this leads to enhanced anisotropy in photons. 

The above results hold for both orders of flow. It is seen that the decay 
dilutes the elliptic flow more than it dilutes directed flow. 

3.2 Estimation of tt° flow 

While the effect of decay is seen to depend on the values of multiplicity and the 
amount of initial flow, the ratio of anisotropy observed in photon distributions, 
v m(l), to that in the parent tt° distribution, v™(ir ), is seen to scale with 
the parameter Xm- Fig. [3] shows this ratio as a function of Xm for a large 
combination of values of anisotropy and multiplicity. 




1 10 1 10 



Figure 3: Ratio of estimated value of anisotropy in photons to the initial value 
in 7T° as a function of measured anisotropy parameter \i f° r different values 
of anisotropy and multiplicity (a) for m — 1, directed flow and (b) for m = 
2, elliptic flow. Symbols denote different initial tt° anisotropics. Logarithmic 
x-axis is used for better visualisation of the results at lower Xm values. The 
error bars are statistical. 
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We have parametrised the observed scaling behaviour by a relation : 

^(7) a 



v in (n°) (x - b)' 



+ c (3) 



The relation was motivated by the fact that the effect of two particle corre- 
lations scales as 1/M [[TJ]] and Xm is proportional to \[M. If the exponent is 
chosen a free parameter, the best-fit value varies between 1.9 and 2.1, depend- 
ing upon the number of points used in the fit. 

The continuous curves in Fig. [5] show the fitted functions for both directed 
and elliptic flow. The values of the constants for different values of m are 
given in Table [1|. The values of b reflect the lower limit of x upto which the 
anisotropy in 7r° can be deduced by this method. The lower value of b for m=2 
implies that elliptic flow in neutral pions can be studied for lower anisotropy 
and multiplicity, c is the limiting value of the fractional anisotropy remaining 
in photons. 



Table 1: Fitted values of the constants of the scaling relation. 



Order (m) 


a 


b 


c 


1 


0.075 


0.473 


0.801 


2 


0.031 


0.286 


0.619 



It is noteworthy that a simple parameterisation can be used up to quite 
low values of x- The parameter Xm depends on anisotropy and multiplicity 
(Xm = v m \/2M) and can be determined from experimental data using eqn. (0). 
Combining this with the anisotropy values v m determined from the experimen- 
tal data on photon distributions, one can deduce the anisotropy present in the 
parent 7r° distribution using the above scaling relation. 

To check the efficacy of the method, we simulated another data set with 
v tn = 0-02 and = 0.04 and for a range of multiplicities varying between 250 
and 3000. This data set was analysed as earlier and the value of anisotropy in 
neutral pions was estimated using the scaling relation. The estimated values 
are close to the input values, the ratio of the two deviating from unity only for 
the cases when both the multiplicity and the initial flow are small. 

The error on the deduced 7r° anisotropy will be small for large values of x-, 
as is evident from the local slope of the curve in Fig. ^. The percentage errors 
on the deduced anisotropy values are shown as a function of multiplicity in 
Fig. [| for both orders of flow. These have been calculated by taking the statis- 
tical errors in f m (7) and Xmi both experimentally measurable quantities. The 
errors are generally within 10%, except for elliptic flow at low multiplicities. 
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Figure 4: The percentage error on the value of 7r° anisotropy, deduced using 
the scaling relation, as a function of multiplicity of photons for (a) directed 
flow and (b) elliptic flow. The circles represent initial anisotropy v m (7r°)= 0.02 
and squares represent v m {ji°)= 0.04. 

4 Summary 

A set of events having built-in azimuthal anisotropies has been analysed using 
techniques employed for the analysis of experimental data. The anisotropies 
are introduced in both charged and neutral pion distributions. Photon dis- 
tributions are obtained from the 7T° distributions. The pion distributions are 
generated to have the desired multiplicities in the range 100 to 4000 in one 
unit of pseudo-rapidity. 

Analysis of charged particle distribution leads to estimates of anisotropy 
values in close agreement with the initial anisotropy present in the event. 
Due to the effect of decay, the values obtained from photon distributions are 
different, being lower, or higher, than the initial anisotropy depending on the 
value of the initial anisotropy and multiplicity. 

The ratio of estimated anisotropy in photon distributions to the initial 
anisotropy in n distributions is seen to scale with a parameter x which can 
be determined from experimental data and depends on both anisotropy and 
multiplicity. This scaling behaviour provides us with a simple relation to de- 
duce the initial anisotropy in 7r° distributions by the experimental study of 
photon distributions. The anisotropy is obtained with varying precision which 
improves when either flow or multiplicity increases. 
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